Simple model of self-supported deformed states of isolated atoms 
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We propose a simple three-body model of an atom in which one electron on a circular Rydberg 
orbit is treated as an independent particle and the remaining core electrons are collectively described 
as a single object. Within this model we predict the existence of stable deformed states of atoms. 
The deformation is generated by a bootstrap mechanism. The atomic core is polarized by the excited 
electron and the induced dipole moment keeps this electron localized. The deformed stable states of 
the atom are similar to the Trojan states observed in recent experiments. However, in the present 
case the breaking of the rotational symmetry does not require the presence of external fields. 



I. INTRODUCTION 

The prediction in Ref. [l[ of nonspreading wave pack- 
ets describing Rydberg electrons driven by a circularly 
polarized wave was recently fully confirmed in experi- 
ments 0, 0] ■ The dynamics of these nonstationary but 
stable electronic states resembles the dynamics of Trojan 
asteroids in the Sun- Jupiter system. The Trojan states 
of atoms are deformed due to the action of the driving 
field. In this article we propose a simple model to study 
a possible existence of deformed states of multielectron 
atoms not subjected to external forces. The deformation 
is caused by the interaction of a Rydberg electron with 
collectively described core electrons. The interaction of 
the Rydberg electron with the core has been studied by 
Clark and Greene Q who described this interaction by a 
phenomenological coupling of the angular momentum of 
the electron to the angular momentum of the ion core. 
They were able to account successfully for subtle de- 
tails of spectroscopic data but the deformation of the 
atoms has not been the subject of their interest. In our 
model we keep the pure Coulombic coupling of the Ryd- 
berg electron to the core but we describe the interaction 
of the core with nucleus by harmonic forces. Recently 
Kalinski et al. [H] have shown in the simplest case of 
the helium atom that due to the core polarizability de- 
formed localized states can exist in isolated atoms. The 
mechanism responsible for the deformation was clearly 
identified in Ref. @ , where we can find the following de- 
scription: "The oscillations of the inner charge generates 
a stabilizing field for the outer electron". The analysis 
presented in Ref. [1] was specifically designed for helium 
and its extension to many-electron atoms would be very 
complicated. We propose a simple model that captures 
the essential features of such special states of atoms. In 
our model one highly excited electron acting on all inner 
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electrons creates a dipole moment and the induced dipole 
localizes the outer electron on its orbit. 

Our aim is to give an intuitive picture of deformed 
atomic states that explains the nature of these states and 
could enable one to pinpoint their characteristics with- 
out the use of more elaborate methods. The purpose of 
the present article is, therefore, to show that the mech- 
anism proposed in Ref. [5j] can produce deformed states 
of many-electron atoms. There is, however, a crucial dif- 
ference between our work and that of Ref. In our 
case the mutual interaction of the core electrons and the 
outer electron leads to a stable configuration. The stabil- 
ity analysis in our model can be carried out analytically. 
The correlations of the electron and the core coordinates 
can be clearly exhibited. The strong point of our analysis 
is that many steps can be performed analytically. 

There is an essential difference between our stable 
electronic states and the ones supported by an external 
field. Since the dynamics of an isolated atom is governed 
by a rotationally symmetric Hamiltonian, the deformed 
states can only occur due to broken rotational symme- 
try. This mechanism is well known in nuclear physics 
where the breaking of rotational symmetry leads to de- 
formed states of nuclei ^, 0| . The analogy between the 
deformed states of nuclei and possible deformed states of 
atoms was pointed out by Iwai Q. 



II. THE MODEL 

To describe intricate dynamics that leads to a de- 
formed atom, we introduce a greatly simplified model of 
a multielectron atom. Within this model we can identify 
atomic states in which the repulsion between one excited 
electron and the core causes the deformation. This lonely 
electron will be described as usual but the remaining core 
electrons occupying orbits close to the nucleus will be 
replaced by a single effective particle endowed with the 
mass and the charge of all the core electrons. The inter- 
action of this effective particle with the far away electron 
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will be kept Coulombic while the interaction with nu- 
cleus will be modeled by a harmonic force with a spring 
constant reflecting the polarizability of the core. The de- 
scription of the strongly bound electrons in terms of a 
harmonic force and the corresponding oscillator strength 
parameters has been successfully used in the early days 
of quantum mechanics. Here we push this description 
further by applying the combined oscillator model to all 
core electrons. Our model cannot give a precise descrip- 
tion of the deformed atom but owing to its simplicity 
gives a clear picture of the mechanism that is responsi- 
ble for the symmetry breaking and for the stability of 
deformed states. 

The Hamiltonian in our model of the atom has the 
following form: 
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where — /Attcq. The first part of this Hamiltonian 
describes the kinetic energies of the nucleus and the sin- 
gled out electron, and also their mutual Coulomb interac- 
tion. The second part describes the kinetic energy of the 
atomic core, its harmonic interaction with the nucleus, 
and the Coulomb interaction of the core with the elec- 
tron. Our model may seem to be best suited to describe 
the alkali atoms for which the distinction between the 
core and the valence electron is well defined. However, 
the polarizability of alkali ions is very small. Therefore, 
as an example, we have chosen magnesium, whose ion 
has the same configuration of electrons as sodium. 

Translational invariance allows for the separation of 
the center of mass motion. This is achieved by the fol- 
lowing canonical transformation: 
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where M — rUe + nic + m„. The inverse transformation 
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The Hamiltonian ([1} expressed in terms of the new vari- 
ables split into the sum H = Hem + ffrci. Choosing the 
coordinate system in which the center of mass is at rest, 
we may drop iJcm = and keep only the Hamil- 

tonian of the relative motion -ffid. 
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where /ig and /ic are the reduced masses of the electron- 
nucleus and the core-nucleus systems. 
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The system described by the Hamiltonian H^-^i has a 
stable stationary configuration. This state in the rotat- 
ing frame becomes an easily identifiable state of static 
equilibrium. In the frame rotating with the angular ve- 
locity vector ijj the Hamiltonian acquires an additional 
term — a;-(r x p + Rx P) describing the inertial forces. 

From now on we shall express all quantities in atomic 
units. The Hamiltonian in the rotating frame H has the 
form. 



p2 p2 p.p ^ 



2fj.e 2fic mn \r\ \r - R\ 
k 

-R^ - u}-{r X p + Rx P). 



(6) 



The canonical equations of motion generated by the 
Hamiltonian H are: 
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For a given atom all masses, Z, and k are fixed and the 
dynamics of the system is controlled by just one param- 
eter oj. This system has two constants of motion: the 
total energy and the total angular momentum. 
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III. EQUILIBRIUM CONFIGURATION 

The system is in equilibrium when all time derivatives 
vanish — all bodies are at rest. The mere existence of an 
equilibrium configuration in rotating frames is not a spe- 
cial property of our model. Such configurations exist for 
all atoms with pure Coulombic forces. They correspond 
to the regular, equidistant positions of electrons. How- 
ever, these equilibria are highly unstable. In our model 
the stability is achieved by a collective description of the 
core electrons. 

In the equilibrium configuration described by 
Eqs. d7a|) - (j7d|l the excited valence electron, when viewed 
from the laboratory frame, moves on a circular Kepler 
orbit with frequency u! around the center of mass and 
the core particle and the nucleus move on much smaller 
circular orbits. After the elimination of momenta, we 
obtain the following two equations for the equilibrium 
position vectors: 
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where ± denotes the component perpendicular to a;. 
These equilibrium conditions express simply the balance 
of forces acting on each particle. It follows from these 
equations that Tq and Rq are aligned and lie in the plane 
perpendicular to u). Since the core is repelled by the 
electron while the nucleus is attracted, these two vectors 
point in opposite directions. We shall choose the x axis 
in the direction of Tq and denote by xi and —X2 the x 
components of Tq and Rq. These parameters represent 
the distances from the electron to the nucleus and from 
the core to the nucleus. Eqs. (|5al) - ((8b| rewritten in terms 
of a;i and X2 form a set of two coupled third-order equa- 
tions, 

Z ^'^e-rrinXi - mcmnX2 . 

— = kx2+u; — , (9a) 
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The first equation is obtained as a difference of equations 
(Hal) and ([gb| . 

Let us denote by 5 the relative distance from the core 
to the nucleus, 



6 — X2/X1. 



(10) 



After dividing Eq. (l9a)) by Eq. (|9b)) we obtain 
Z{l-\-5y' kMS + Lo'^irriemn — mcrrinS) 



Z — I kM6 — uP' [merric + {rricmn + memc)S] ' 

(11) 

and this leads to a single third-order equation for S. Full 
solution of the equilibrium conditions can be obtained 
in the general case. Here we present only the results 
obtained in the limit of an infinite nucleus mass. We 
could proceed without making this approximation, but 
in the case of a many electron atom the results would 
differ by a small fraction of a percent. In this limit the 
equations for xi, X2, and 6 read 



— =u^xi + [k- {Z -l)uj^]x2, (12) 
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where we explicitly used the fact that the core mass mea- 
sured in units of the electron mass is Z — \. It follows 
from Eq. that the value of 5 tends to zero together 
with uj, 



Z - 1 



(15) 



Since 5 must be positive, Eq. (jl4[) imposes an upper 
bound on w, namely {Z — l)cj^ < k. This bound simply 
means that when the rotation is too fast the centrifugal 
force will overcome the spring tension and the core will 
fly away. On physical grounds this upper bound is a bit 
too large. When the distance from the core to the nucleus 
becomes equal to the distance from the electron to the 
nucleus, our model becomes clearly inapplicable. This, 
according to Eq. ([T2|) . leads to a slightly more stringent 
condition 



3Z + 1 



(Z- l)(3Z + 2) 
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Equation (|14p can be easily solved for w. 
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Now, dividing both sides of Eq. (|12p by xi we obtain the 
distance Xi from the electron to the nucleus as a function 



of cli, fc, and Z, 



<:i{uj,k,Z) = 



-,1/3 



or as a function of (5, fc, and Z , 

■(Z-l)(l + 5 + 2Z(52 + Z^3)ii/3 



fc, Z) 



5{l + 5Yk 



(20) 
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The distance X2 from the core to the nucleus is obtained 
by multiplying xi by 5. When w — > we obtain 



Xi{uj, k, Z) 



1 



,2/3- 



(22) 



This result is in agreement with the classical formula for 
the radius of a high hydrogenic Rydberg orbit because 
for such orbits the only effect that the core has on the 
electron is the screening of the nuclear charge. In the 
equilibrium configuration, the z-component total angular 
momentum along u) is 



: to [xi{uj, k, Z)"^ + {Z - l)x2(w, fc, Zy] 



(23) 



IV. 



PHYSICAL INTERPRETATION OF THE 
EQUILIBRIUM CONFIGURATION 



The results of the previous section confirm our intuitive 
picture of the deformed atom. As might be expected, in 
the equilibrium configuration all three objects, the core, 
the nucleus, and the electron, lie on a straight line. The 
core and the electron are on the opposite sides of the nu- 
cleus. In planetary terminology, the equilibrium config- 
uration in our model corresponds to the Lagrange point 
L3 — all bodies are aligned and the planets are on the op- 
posite sides of the Sun. Unlike its planetary counterpart, 
however, our equilibrium configuration is stable (for suf- 
ficiently small a;). In the next section we prove the sta- 
bility of motion near the equilibrium by calculating the 
characteristic frequencies in the quadratic approximation 
to the Hamiltonian. 

All calculations in this article will be done for the mag- 
nesium atom: rric = llrrie, vrin = 43710 mg, and Z = 12. 
To complete this set of data we will need the value of the 
spring constant k. This can be determined from the po- 
larizability. We can relate the parameter k in our model 
to the polarizability a, defined in the standard manner 
as the ratio of the induced electric dipole to the applied 
electric field. The induced dipole is {Z — l)eR and the 
electric field E of the electron that creates this dipole by 
displacing the core away from the nucleus is 
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FIG. 1. The relative distance 5 from the core to the nucleus 
as a function of uj for magnesium. 




FIG. 2. The distances between the electron and the nucleus 
x\ and between the core and the nucleus X2 as a function of 
uj for magnesium. Distances are measured in atomic units. 



Using the values at the equilibrium, we arrive at the for- 
mula 
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The right-hand side is a function of the model parameters 
and Lo. In order to obtain the static polarizability we 
should take the limit when w — or 5 — 0. Therefore, 
we may use the asymptotic formula. 



1 



k5 



(26) 



obtained from Eq. ([2T|) when d ^ Q. In this limit we 



5 



arrive at an approximate relation between the static po- 
larizability and fc, 



Aneok 



(27) 



This formula can also be obtained directly from the 
following observation. In the limit when |ii|/|r| — 7> 
and m„ — > oo the Hamiltonian (|4]) describes two non- 
interacting systems: an electron in the Coulomb field 
and an oscillator with mass {Z — l)me and the spring 
constant k. The displacement of an oscillator of charge 
{Z — l)e under the action of the applied electric field E 
is X ^ {Z — l)eE/k. Thus, the polarizability defined as 
the ratio of the induced electric dipole {Z — l)ex to the 
electric field [Eq. ([Ml)] is given by the formula Eq. (P7|) . 

Using the value of the rationalized static polarizability 
a of the magnesium ion [9] , 



(5 = 34.62 a.u. = 5.13 x 10 
we obtain the following estimate 

k w 3.5 a.u. = 5442 kgs" 



(28) 



(29) 



The characteristic energy of the core oscillations corre- 
sponding to this value is h^/kjrric = 15.3 eV, a fairly 
reasonable number. Having fixed the parameter k of our 
model, we can proceed with detailed calculations. 

All numerical calculations and figures in this article 
were done with the use of MATHEMATICA [13]. The 
function d{uj, Z) for magnesium is depicted in Fig.[TJ The 
value of 5 — the relative distance from the core to the 
nucleus — is a good measure of the atomic deformation. 
In Fig. [2] we plotted the values of the equilibrium coordi- 
nates Xi and X2 expressed in Bohr radii qq. At the value 
[Eq. ([T6| ] of oj they become equal. 



V. STABILITY OF A DEFORMED STATE 

We will investigate the stability of the equilibrium 
state in our model by linearizing the equations of mo- 
tion near the equilibrium configuration. This standard 
method worked very well in the studies of Trojan states 
[J, [nl, [ll] . We shall begin with the linearized equations 
of motion obtained from Eqs. ([7|) in the limit of infinite 
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where the two-dimensional vectors cci, /3i, 0:2, /32 de- 
scribe small deviations from the equilibrium positions 
and momenta of both particles in the xy plane, while 
Ci, 0-1X2, ^2 describe the deviations along the z axis. The 
unit vector parallel to the line passing through the equi- 
librium positions (in our case the x axis) is denoted by 
m. Note that in this approximation the motion in the z 
direction decouples from the motion in the xy plane. 

The right-hand side of the linearized equations of mo- 
tion can be encoded in the following two matrices: 
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FIG. 3. Frequencies of the small oscillations in the xy plane 
as functions of the rotational frequency lo for magnesium. 
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FIG. 4. Frequencies of the small oscillations in the z direction 
as functions of the rotational frequency uj for magnesium. 



The matrix Mxy acts on the eight-dimensional vec- 
tor {a.1, /3i,a2, P2}, while the matrix Aiz acts on the 
four-dimensional vector {Ci, cri, C2j o'2}- The eigenvalues 
of these matrices determine the frequencies of normal 
modes. 

The matrix A4xy gives four pairs of frequencies — the 
members of each pair differ only in sign. One pair consists 
of zeroes — a result of the broken rotational symmetry. 
By analogy with field theory, the motion corresponding 
to the zero frequency may be called the Goldstone mode. 
The positive partners of the remaining frequencies are 
plotted in Fig.[3]as functions of the rotational frequency uj 
in atomic units. The smallest eigenfrequency wi vanishes 





FIG. 5. Oscillations of the core (dashed line) and the elec- 
tron in the xy plane around their equilibrium positions for 
magnesium at cj = 0.06 a.u. corresponding to the slow mode, 
uJi = 0.041 a.u. The size of the core orbit was increased ten- 
fold to make it visible. In the slow mode the core and the 
electron move in the same directions along the x axis but in 
the opposite directions along the y axis. 



at uj =0.1026 a.u. and for larger values of ui it becomes 
complex. This sets the stability limit of the equilibrium 
state. The maximal value of S is 0.022. In the stability 
region, the value of the classical angular momentum 
changes from 2. 8/1 to 00. 

The matrix A4z gives two pairs of frequencies. The two 
positive frequencies are plotted in Fig. |4l In the limit, 
when a; — > the nonvanishing frequencies uj2, ^3 and W5 
take on the value y^k/me{Z — 1). This frequency can be 
viewed as the characteristic frequency of core oscillations. 
For magnesium this value is 0.5637 a.u. 



VI. CLASSICAL ORBITS 

We begin with the description of small oscillations cor- 
responding to different eigenmodes. To a good approx- 
imation, these small oscillations are harmonic — the or- 
bits are elliptical. In Figs. [5] and [6] we plotted small 
oscillations in the xy plane of the core and of the elec- 
tron around their equilibrium positions for lo = 0.048 
a.u. They were obtained by calculating numerically the 
three eigenvectors of the matrix Mxy corresponding to 
the eigenfrequencies wi, W2, and 0^3. In all three modes 
the motions of the electron and the core are strongly cor- 
related. In the slow mode (Fig. [S]), corresponding to Wi, 
the amplitude of the electron oscillations is much larger 
than the amplitude of the core oscillations while for the 
fast modes (Fig. [6]) the amplitudes of the core oscillations 
dominate. 

Classical orbits of the core and the electron without 
the approximation of small oscillations can be easily ob- 
tained from the numerical solutions of the full equations 
of motion [Eqs. ((7a)) - (j7dl) ]. These solutions show the ro- 
bustness of our system far beyond the regime of small 
oscillations. For a random choice of the initial conditions 
all the modes become excited, the core oscillates rapidly 
(see Fig. [7]) but the electron follows similar smooth tra- 
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FIG. 6. Oscillations of the core (dashed line) and the electron 
around their equilibrium positions for magnesium at a; = 0.06 
a.u. corresponding to the two fast modes, 102 = 0.506 and 
Lo-i — 0.626 in a.u. for the same parameters as in Fig. (5] In 
all figures the norm of the amplitude is the same as in Fig. [S] 
In both modes the core and the electron move in opposite 
directions along the x axis. The modes differ in the motion 
along the y axis. In order to keep the same overall size as in 
Fig. [5] (assuming the same norm of the amplitude), the overall 
size was scaled up fivefold. 
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FIG. 7. The trajectories in the xy plane of the core (left) and 
the electron (right) in the regime of small oscillations when 
all modes are excited. The parameter oj = 0.048 corresponds 
to total angular momentum 3ft. At the start the electron was 
displaced to the right from the equilibrium position (marked 
by the dot) by 0.2 a.u., while the core started at the equi- 
librium position. The time lapse is equal to the period of 
electron as it covers the elliptic orbit. 



jectories as in the regime of small oscillations. When the 
departure from equilibrium is large, the electron trajec- 
tory is still regular and stable but it acquires a banana- 
like shape (see Fig. [S]). 



VII. QUANTUM MECHANICS OF A 
DEFORMED ATOM 

In the description of a deformed atom in quantum me- 
chanics we shall start from the same model Hamiltonian 
in the rotating frame. Unfortunately, we cannot fol- 
low the same path as in the description of the Trojan 
states [l|. 111] and use the approximation of small oscil- 
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FIG. 8. The trajectories in the xy plane of the core (left) 
and the electron (right) during one cycle of the periodic elec- 
tron motion, outside of the regime of small oscillations. The 
banana-shaped orbit of the electron shows that the system is 
fairly robust. The electron was displaced from its equilibrium 
position (marked by a dot) by as much as 5 a.u. in the hor- 
izontal direction and was given the initial velocity 0.12 a.u. 
in the direction of rotation. The core was initially displaced 
from its equilibrium position by 0.1 a.u. The total angular 
momentum of the system is 3h. The nucleus is in the center 
of the coordinate system. 



lations. This is due to the fact that the characteristic 
frequencies of the electron motion cji and 0:4, as shown 
in Figs. [3] and H] are tiny. Hence the extension of the 
electronic Gaussian wave function will exceed the range 
of small oscillations. Therefore, to find the quantum sig- 
nature of the lateral localization of the electron vis-a-vis 
the position of the core we must use a different approach. 
We shall approximate the interaction potential between 
the electron and the core by the first two terms of the 
expansion in the ratio |ii|/|r|. 



Z - 1 
\r-R\ 



Z -1 



+ (Z-1 



Rr 



(34) 



This approximation is justified, since we aim at the study 
of electron states localized far from the core and it is also 
well supported by the classical theory. In this approxi- 
mation, the Hamiltonian reads: 

H = -W-^+^u:■{rxVr) + {Z-l)^ 



1 



2{Z -1 



-Ar + -R^ +iuj-{R X Vi 



(35) 



We took the limit of an infinitely heavy nucleus and we 
used atomic units. Even after these simplifications, an 
exact solution of the quantum problem is a hopeless task 
because the variables cannot be separated. However, we 
can simplify the problem following the hints provided by 
the classical theory. 

The part of this Hamiltonian involving all the terms 
describing the core has easily calculable eigenstates. The 
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FIG. 9. Plots of the electron density along the x axis (left) and 
of the density in the xy plane (right) for the superposition of 
n = 5 and n = 4 states. The dot indicates the position of the 
nucleus. The interpretation of contour lines is made possible 
with the plot of the effective density p along the x axis in 
the equatorial plane. The coordinates x and y are measured 
in atomic units 



lowest energy state is described by a shifted Gaussian 
(i){R,r) = iVexp 



'^-{R + gf +i{Z-l)u^-{Rxg) 



where 



7= V^(Z-l), 

1 ( X V z 

9= ' ^ 



k 



(Z- 1)^2 



(36) 

(37) 
(38) 

(39) 



and N = {'^/ttY^'^ is the normalization factor. As in the 
classical case, we have chosen the z axis along oj. The 
core wave function still depends on the electron coor- 
dinates through the vector g. Since the motion of the 
electron is slow as compared with the fast oscillations of 
the core, as shown in Fig. [71 we may approximate the 
wave function \E'(i2, r) of the whole system by a Born- 
Oppenheimer-like product. 



^{R,r) 



{R,r)^{r) 



(40) 



Within this approximation we can confirm our classical 
picture of the deformed atom in which the electron and 
the core are localized on the opposite sides of the nucleus 
since the expectation value of the scalar product il-r is 
strictly negative, 

j (fR(fr\(j){R, r)^{r)fR-r = - j (fr\ij{r)\'^g-r. 

(41) 

We carried out here the Gaussian integration over i?, 
after a change of variables i2 — >■ R — g. 



FIG. 10. The same plots as in Fig. [9]but for the superposition 
of n = 6 and n — 5 states. 



Next, we will determine the shape of the electron wave 
function while the position of the core R is kept fixed at 
its classical equilibrium value R = (— a;2,0, 0). This as- 
sumption is based on the classical analysis which showed 
that the core undergoes rapid oscillations in the close 
vicinity of the equilibrium position while the outer elec- 
tron moves much more slowly (see Fig. [7]). Under these 
assumptions the Schrodinger equation for the electronic 
wave function "tp^r) is 



where 



Ho 

Hi 



{Z-l 



1 

Rr 



iu)-{r X V), 



(42) 

(43) 
(44) 



In the region of large r, where we expect electron lo- 
calization, the dipole term Hi is very small and can be 
treated perturbatively. The eigenfunctions 'ipnim{'r,9,4>) 
of Hq are well known. The only possibility for a small 
perturbation to produce a large modification of the wave 
functions ^{r) is the occurrence of a resonance. Due to 
the presence of the rotational term —ujMz in Hq, there is 
a possibility of such a resonance. This is a typical case of 
level crossing and it happens when the rotational contri- 
bution ujMz to the Hamiltonian matches the difference 
between the hydrogenic energy levels En = —Xjlr? in 
the rotating frame. 



(45) 



In such a case we have to deal with the degenerate per- 
turbation theory. The Hamiltonian El\ has nonvanishing 
matrix elements between pairs of states that satisfy the 
dipole selection rules AZ = ±1, Am = ±1. For definite- 
ness, let us choose n\ > n2- Without loss of generality, 
we may assume that Am ~ mi — m2 = 1 and we obtain 
the resonance value of the rotational frequency. 



n\n2 — ^Tii 



E„ 



(46) 
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The opposite sign of Am would only give the reversed 
direction of rotation. According to our classical analysis, 
the deformation will be the largest when the electron 
stays near the equatorial plane and is sufficiently far from 
the nucleus. This means that I and m should be the 
largest possible. This observation goes very well with the 
fact that matrix elements of Hi between the hydrogenic 
wave functions ipnim{'>'id,4') ^re maximal for the largest 
values of I and m. Last but not least, we must stay 
within the limits of our version of the Born-Oppenheimer 
approximation and of the dipole approximation — both 
require the outer electron to be far from the core. On 
the other hand, the interaction between the core and the 
electrons weakens with the increasing distance between 
them. 

We shall build the quantum state of the deformed atom 
from Coulomb states with maximal values of the quan- 
tum numbers I and m, for a given n. Under this as- 
sumption, the subspace of states with the same eigen- 
value of the unperturbed Hamiltonian, for the rotational 
frequency w satisfying Eq. (|46p. is two-dimensional. For 
each n, this subspace is built on two hydrogenic states: 
\n, n — 1, n — 1) and \n — 1, n — 2, n — 2) . The degen- 
eracy is removed by the interaction term represented by 
the matrix 



dn 
dn 



(47) 



where d„ is the matrix element of the interaction Hamil- 
tonian Hi taken between the hydrogenic wave functions, 

R r 

dn = {Z — l){n,n — l,n — 1\- — —In — l,n — 2,n — 2) 

X 

= —X2(Z — l){n, n— l,n— l\-—-^\n — — 2,n — 2). 

(48) 

The normalized eigenvectors of the matrix (|T7)) have the 
form: 

|n±) = (|n,n - 1, n - 1) ± |n - 1, n - 2, n - 2)) . 

(49) 

These states correspond to the values of the quasienergy 
En - {n- l)u} ± dn- 

Superpositions of states with different values of the 
magnetic quantum number always lead to wave packets 
of the electron density on the orbit. However, in general, 
such wave packets will not be stable. In our model, due to 
the interaction with deformed core, the energy difference 
stabilizes the deformed states. 

The superposition [Eq. (|49| ] with the plus sign, gives 
the electron density localized on the same side of the nu- 
cleus relative to the average position of the core described 
by R. This state has a higher quasienergy {dn > 0). The 
superposition with the minus sign leads to the localiza- 
tion of the electron on the opposite side of the nucleus. 
Therefore, it corresponds to the lower quasienergy. This 



quantum state, together with Gaussian state of the core 
(j>{R,r), is a counterpart of classical localized orbits of 
the electron and the core — it gives a quantum descrip- 
tion of a deformed atom. 

To illustrate these phenomena we choose the follow- 
ing superpositions, including their time-dependent phase 
factors 

|5-) = exp[-z(i?5 - 4t^54 " d5)i]^ (|5, 4, 4) - |4, 3, 3)) , 

(50a) 

,4)). 
(50b) 



|6-) = exp[-t{Ee - 5^65 - ^e)^]^ (|6, 5, 5) - |5, 4, 4)) 



The electron densities are depicted in Figs. |9] and [TO] 
We see that the localization weakens with the increase of 
the principal quantum number. This was to be expected 
since a larger n means a weaker interaction with the core. 

In what follows, we shall focus on the state |5— ). The 
matrix element in Eq. (|^5)) for n = 5 is equal to —0.001. 
According to (j46|). in this case the value of a; = UJ54 is 
1/32 — 1/50 and X2 = 0.0079. Therefore, the quasienergy 
for this state is decreased by (^5 = 9 x 10~^. 

It is instructive to transform the stationary state in the 
rotating frame to the (inertial) laboratory frame. This 
transformation is accomplished simply by multiplying ev- 
ery eigenstate of belonging to the eigenvalue hm by 
the phase factor exp(— iwmt). This means that the state 
vector (|50ap transforms into 

|5-) = e-*(^s-4c.-d5)t 1 (e-*4<^*|5,4,4) - e-*3'^*|4, 3, 3)) 

- e-'^^^-"^^)*-^ (|5, 4, 4) - e-*|4, 3, 3)) . (51) 

The electron probability density in this state is given by 
the formula: 

p{r,9,cb,t) = i (|V'5,4,4(r,^)P + |V'4,3,3(r,^)P) 

- V'5,4,4(''> &)'4'4,3,3{r, 6) cos{(p - UJt), (52) 

where we denoted by ipnim (f, 0) the hydrogenic wave 
function without the phase factors exp(irniy9) and we have 
taken into account the fact that these functions for the 
maximal values of I and m are real. Thus, as expected, 
the superposition (|50ap transformed to the laboratory 
frame describes an electron wave packet revolving around 
the nucleus with the frequency (I46p . 



VIII. RADIATIVE DECAYS AND THE 
ROTATIONAL FREQUENCY SHIFT 

The quantum state |5— ) of the deformed atom has only 
two decay channels with the emission of dipole radiation. 
The selection rules and the quasienergy conservation al- 
low the transitions to the states |5-|-) and |3,2,2). In 
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the first case, the radiated photon has the quasienergy 
— 2(i5 and it is linearly polarized along the y direction. 
The dipole matrix element (5 + |j/|5— ) is equal to — 8.35i. 
In the second case, the quasienergy of the photon is 
2E4^ — — — and the radiated photon is emit- 
ted with circular polarization. The dipole matrix element 
(3, 2, 2\{x — iy\^—) is equal to a/24.74. These results must 
be now transformed to the laboratory frame. 

Deformed states of the atom offer a unique opportunity 
to study the effects of the rotational frequency shift pre- 
dicted some time ago [l^. These effects up to now have 
been seen only in molecules (l^ . Rotational frequency 
shift is an analog of the Doppler shift. To account for 
the Doppler shift, when comparing the energies of a pho- 
ton with momentum hk measured in two reference frames 
moving with respect to each other with velocity v, one 
has to add the scalar product v k. To account for the 
rotational frequency shift, when comparing the energies 
of a photon with angular momentum M measured in two 
reference frames rotating with respect to each other with 
the angular velocity u;, one has to add the scalar prod- 
uct U3 ■ M of the angular velocity and photon angular 
momentum. The rotational frequency shift does not lead 
to any observable phenomena for a state that is rotation- 
ally symmetric since there is no way of telling whether a 
quantum system in such a state is rotating or not. 

The transformation to the laboratory frame does not 
affect the transition matrix elements but the energy of 
the emitted photons acquires an additional term uj = 
E^ — E4. Thus, the spontaneous decay of our state |5— ) 
leads to the emission of a linearly polarized photon with 
the energy E^ — E4 — 2d^ or a circularly polarized photon 
with the energy E^ ^ E3 — d^. The transition rates are 
practically the same as for standard Rydberg states with 
maximal values of I and m. 



IX. CONCLUSIONS 

With the use of a simple model of an isolated multi- 
electron atom we predicted the existence of stable states 
describing a deformed atom. In such states, one excited 
electron is localized as a wave packet on a circular Ry- 
dberg orbit. The core made of the remaining electrons 
is polarized due to the interaction between the electron 
and the core. The deformation requires breaking of the 
rotational symmetry. In classical theory, we were able to 
give exact description of the atomic deformation and we 
proved the linear stability of deformed states. In quan- 
tum theory, we had to resort to a perturbative treatment 
of the interaction between the core and the electron. Lo- 
calized states of the excited electron are special superpo- 
sitions of two circular Rydberg states. In the construc- 
tion of such stable localized states, the essential role was 
played by the phenomenon of level crossing in the frame 
rotating with the electron. 

We cannot give a prescription on how to produce 
the special superpositions describing a localized electron. 
However, the calculated redshifts in spontaneous decays 
of those states could be used as their signature. Unfortu- 
nately, our model is so crude that its predictions cannot 
be precise but we hope that they are of the right order 
of magnitude. 
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